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I.  SUMMARY 


During  the  course  of  this  research,  three  technical  areas  have 
been  addressed: 

(i)  eddy-meanf low  energetics 

(ii)  model  extension  to  irregular  geometries,  and 

(iii)  flows  around  seamounts  * 

In  the  first  area,  three  models  have  been  developed  to  study  eddy-mean- 
flow  energetics  in  ocean  circulation.  The  formulation  of  these  models 
is  unique  and  it  is  felt  that  they  are  the  first  truly  applicable  to 
regional  energetics. 

Similarly,  NORDA  modeling  capabilities  have  been  considerably 
extended  by  permitting  the  modelling  of  irregular  ocean  basins.  This 
extension  from  simple  rectangular  geometries  has  greatly  improved  the 
realism  of  the  model  and  its  behavior.  - 

The  seamount  studies  are  specifically  directed  at  an  examination 
of  measured  circulation  perturbations  near  a  seamount  topography  and  a 
comparison  with  predicted  patterns.  These  comparisons  provide  data  on 
model  ability  to  predict  flow  fields  and  also  on  the  accuracy  of  pre¬ 
dictions. 

Each  of  these  areas  is  discussed  in  the  remainder  of  this  report. 
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II.  EDDY-MEANFLOW  ENERGETICS 


Three  software  packages  are  developed  and  tested  for  the  study  o* 
eddy-meanflow  energetics  in  ocean  circulation.  They  are: 

(i)  one-layer  g' -model  energetics 

(ii)  two-layer  free-surface  energetics  (large  core  version) 

( i i i )  two-layer  free-surface  energetics  (small  core  version) 
(ii)  and  (iii)  can  also  be  used  for  rigid-lid  model, 
a  separate  preprocessor  is  provided  for  this  purpose. 

The  ocean  models  when  subjected  to  steady  forcing  functions  usually 
reach  a  steady  or  statistically  steady  state,  provided  some  dissipation 
mechanism  is  present  and  integration  is  carried  out  long  enough.  For 
statistically  steady  cases,  a  meanflow  can  be  defined  as  the  average  of 
the  flow  field  over  a  period  of  time  in  the  final  trendless  state.  The 
eddy  field  is  defined  as  the  difference  between  the  instantaneous  flow 
and  the  meanflow.  Energy  equations  are  then  derived  for  the  meanflow 
and  the  eddy  field;  they  relate  the  time-rate-of-change  of  various  energy 
measures  (Kri),P(I1,Ke,Pe)  to  the  energy  transfers  between  them.  A  schematic 
energy  diagram  is  often  constructed  to  quantitatively  show  the  flow  of 
energy  between  the  different  variables. 

Several  energetics  studies  have  been  done  in  the  research  of 
mesoscale  eddies  and  general  circulation.  Most  of  them,  however,  are 
concerned  only  with  the  global  energy  balance,  or  integral  properties  for 
the  entire  basin.  Our  formulation,  we  believe,  is  the  first  for  layered 
ocean  models  that  has  correct  Reynold's  stress  terms  and  the  capability  of 
regional  energetics.  Detailed  derivations  of  the  energetics  equations 
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have  been  delivered  to  Drs.  J.  D.  Thompson  and  H.  G.  Hurlburt  of  NORDA. 
Because  of  the  length  and  the  onerous  mathematical  formulae,  the  deriva¬ 
tions  are  not  presented  in  this  report.  Rather,  we  concentrate  on  the 
usage  of  these  energy  packages. 


A.  One-layer  Energetics  (q'-model) 
(i)  Run  Preprocessor 


/  J3SJ  LIN»G?M.EV1LPR£(J  ,  i  7  X  ~  4  6  A  «LI  .1  Li  , PP T  =  , D 1 4 ) , C M *34 

/  L  I u  IT  U  A  ft  D  =  4  5  0  t  *'  I  ‘!  ■=  2  . 

/  PD  M  ,UStRCAT/'J77/r  S  c /L  [  1  L  1 

/  FO  FTu9FCGl  ,  r;  A ! :  L)  =  2  G  /  3  3  /  i>  >  F  !Ro  =  P  S  ,  o'<  S I  =4  OP  *  r  C“  •»  9  m  s  ,  C  r- M  =  V  P  S 
/  F 0  F T U  9 F C 0  2  »  3  AND  =  2 0 /  3 6  /  F  ,  p  j R  3= R  $  ,  r  ;<  $ Z  =  4  j 0  0  >  L R  E  C  =  j  7  9  6  »  R  C F  R  =  V 2  S 
/  cO  FTo9cO03  i  j  ViD=20  /35  /i> »  -  Vn'j-PS  n:  X  SZ  =4  uG  '.*»  L  C=39Q  a  ,  i  c  r‘i  =  V2  S 
/  9  F  R  FiLEXD  *ErID  =  0'j6B34,uA31=j 
/  FIT  FT  CDF  001 
/  c  IT  F  T  C-  9  r  0  C  2 
/  FIT  FT09F003 
/  9FRn 

/  FO  FT  1 OFOO 1 »  2  4fti)  =  150/2  3  Z  /22  « F  P. G  =  PS  » SZ  =  4C  3  "• > L2"D=  '2  9  16 ,  R C  F'1 -V  2  $ 
/  ASG  fll,  1Y/tNGlLYk.Ts(JS-  =  5  1 1 
/  L NX  LNKTPT  =  ('J,Y,A»S,F,  j,K) 

LITRAKY  11 
INCLUDE  ’VU.PRO 

/  FXGT  1  T=(T,Z), CPT19E=9.io  ■  - 

E.IDP  r 

R4  0=l(.  .j  0.  ,TC=-lf  SPV  «•_=..  wE  +  lc,  ,,  =  3»N3t>  =  36, 

GE  NO 

/  CATV  M  Y/F  46 8  34 Y3  »  ACft  S  =F  1 1  '  F  C  .  2 

/  ETJ 


The  history  tape  is  #006834,  the  output  is  a  long  file  containing 
three  years  of  data  and  resides  on  disk  "EM6834Y3". 

See  program  listing  for  the  meanings  of  the  input  parameters. 
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(ii)  Run  Energetics  Program 


/  JOf  Lld.G33.fc“ElL'>U'i,17  C4S)u  il!  ,*  L  ’  ,  **  P  T  =  CC  ,  1 ,  )  ,  C  T  =  , 

/  LIMIT  aA?<r)  =  ‘ir',MI1'j  =  3 

/  PD  MY  .UStPCAT/ )77/bSu  'Ll  .IL  1 

/  ASG  01,  MY/EGolLrn3t0S,'  =  S.i  ’ 

f  ASG  OrtJ,*'Y/NORDATEK,US;=S 

/  ASG  FT1DF0C1  ,MY/C.'1683mY3,  )SE  =  S  i- 

/  LNK  L'i,sflPT  =  (N,Y,/»,S,F*  i ,  “  )  »  LSP -.C  E=  2C  ? GI 

LIBRARY  01 

INCLUDE  EuGlLr. 

LIBRARY  3  3 J 

/  CD  c T 1 1 c 0 u 1 »  C  A  ■  n  =  4  / 1  j  /  i 
/  FXOT  3  P  T  =  C  Y y  L  )  ,CPT i N  r - y  u  j  0  3 
CINP'JT 

ISTH  =  1,  IMUH=R0  ,  JSTU=4  ,  JV).1  =  48, 

IFLA3  =  1  ,inAY  =  34  0  ,  IS.<IP  =  1 1  ,  nx'.F.C=  >7  , 

CtND 

/  CATV  MY  /  L  R  ASS  ,  AC  VI  =  c  T  *  i  j  j  i 
/  EOJ 


The  output  from  preprocessor  ("EM6834Y3")  is  used  as  the  input  for 

♦ 

this  program.  The  user  should  examine  the  printout  of  the  preprocessor  and 
determine  the  time  domains  of  the  averaging  process  (i.e.  I  SKI P ,  MAXREC, 
and  IDAY).  See  listing  for  meanings  of  input  parameters. 

i 


* 
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(iii)  Switching  Storage  Device  for  Data  File 


/  JOB  LIN.  GDM.T  APl2DS*;  ,1  7  a  ^  ^  K;  ,L  I  iTLl  ,  T=(C,  0,  l)  ,CA  T=1  o 

/  LIMIT  bANT  =  280,*<I‘J  =  : 

/  PD  MY  iUStP.CAT/D77/J5D/L  INSL1 
/  FD  ^»8&'JD=15C'/ZSC/2:; 

/  A  SG  F ,  M  Y / L  fi 6 8  3 4  Y  3  »  U S i.  -  i i : K 
/  CATV  MY/EM6B34Y3. ACNN=r 
/  EOJ 

- /  / - 

/  JOB  LIN.GOB.DSK2TA(JE,:7X  '4bX  '  ,L~  Ji  LI  ,  *<r  T  =  C  C  ,  :>  C  '  T=i 
/  LIMIT  3  A  N  D  =  5  0  »  I  N  =  1 
/  PD  MY,USfcP.CAT/D77/^5 
/  A  SG  F,  MY/EM6334Y3  ,USu-A  v. 

/  CATV  KY/EM6834T3,  ACM-'.=r  ,">TYP  =  7  \r  ■ 

/  EOJ 


The  preprocessor  creates  data  file  "EM6834Y3"  on  disk,  ready  for 
use  by  the  energetics  program.  However,  disk  changes  for  such  a  long 
file  are  quite  high.  It  is  advisable  that  the  user  move  the  file  to  tape 
after  each  use.  The  above  two  short  jobs  switching  the  storage  device. 
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(iv)  Program  Modification 


/  JPF  LINt&0*iiENi:Rur,17<v45)(;  ,Li.l«Ll  ,T’T  =  CC,r',=!),Ctr-i 
/  L  IM  IT  =  > 

/  PD  :1Y,USCRCAT/D77/D5 J/LINTLI 
/  ASG  S 1 »  MY / E  ?JG  1L  YR  S  » 'J 5 F  =  3.1 1 
/  ASG  PI  ,  MY  /  E  N3 IL  YRO  ,IJS =  S  V’ 

/  CIFER  CAOOME.M  =  33K,SEC  =  24'' 

<<  C DP Y  S1,S2 
<<  C DP Y  01,. 

<<  MERGt  Si, TEMP, UPDATE 

<  SELECT  MASK  <" - 


-10,13 

1  I.NDP2*  JSTPt  I-MX  »J  J  ■  i'Jf'Z,!  lAX"*? 

-43 

I.MAXP2=IM 4X+2 
<  SELECT  PREPR3 
-152,152 

83  3  FUR  MAT  (IX  ,  *"3EC  =  '  ,  j.a,5X,  'KE  =  '  ,1  ?  ?.  12. 3,  '  0  _  =  '  ,  1  ”  E  1 2 .  3) 

/  F  TG  F  T  V  E  R  S  =  NX  ,  F  T  VJP  T  -  (  ’J  ,  X  ,  i ,  Y  ,  V' ,  J  )  »  I  N  =  T  f  ' p  ^ - 

/  FD  eTu6sC01f  ..iA*,'n=4/l  >/c 
f  CIFfcR  C  A0DME.M  =  3DK  ,  SCC~2  4 
<<  SDL  1 1  T  E  v  P  ,  S  2 
<<  SPLIT  SYS. 3.130,32 
«  PRIST  S2 

/  CATV  MY/ENG1LYRS, AC  NS -S 2 
/  CATV  M  Y  /  fc  »|  j  1  L  Y  R  .1 ,  A  C  N  '  -  ' 
t  EOJ 


In  the  event  that  the  user  wants  to  make  permanent  changes  to  the 
energetics  package,  the  above  program  provides  the  CIFER  updating 
procedure.  Changes  should  be  bracketed  between  <<  MERGE  and  /FTN. 


B.  Two- layer  Energetics  (small  core  version) 
(i)  Run  Preprocessor 


t  Ll'i.oV'.SC'  V  ■’  2  t  J  ,  l  7  c  y  ,  LI 

/  LI- IT 

/  po  ryjIjsercat/:  77/js  /Lr.v_: 

/  ASG  01  »  o  Y  / S C  J  J ;  L R  J  i  J i>  r  =  S. I 
/  L  fl<  LN<  1PT  =  (-J,  v  f  A  ,  J  ,F  ,  . 

LIBRARY  ”1 


•  V  - 


•’  ?  T  =  (  . 


)  . 


I  4CLIJ  IE  P  R  f:  P 

f  F  0  r  Tu  9  CC  0  1  »  L  V. .)  =  a,  C  / i  i  l'  /i  I  j  ~  l ”  j  =.-> 3  i  ,  $  2  =  4  i  ,  ! 

/  FIT  F  T  0  A  F  C  A  I  ,  £  F- 1 D  =  J  0  2  I  2  i  ,  L  A  -  L  =  " 

t  FT  4  1  PC  C  1  t  A  A  L'=  4C  /  7  /*,-"■»  j  =  "  S  ,  :  $ Z  =  4  . ;  :c 

/  F  0  FT42cCj1  juI,  iJ=4''/7  /*--  ,-«j  J-.J(  =4  -  • 

/  fx:t  ',pt=cy,z  ),cptij-  =  i  , 

cr,°j- 


R  ■!  3  -  1 1  u  .  >  i 

I=L  =  i  ,L  .\ST=3, 

Cth  J 

/  :atv  r /IJ2  j 22 5  »  \C  I  -  =  .: T4  i  - ;  •  i 
/  CATV  •'  Y  /  L  2  f>  2  2  “  .  5  »  A  v  V'  -  ^  I  4  > ' 
'  (,«“ 

/  FIT  FT?*)  r-C  -I  I  ,  E  "  £ «:  =c  0  ;  2  4  3  ,  •  L  = 

f  FACT  AHT=CY,2T,CPTx'.E^  y.-> , 
criF'iT 

BA^lOuC  .  ,-3C  =  -I  V-6=  •  s, 

IcL  =  2iLA>f>T  =  jj 
CE  V~ 

'  CATV  -r/U.'242F.;i,  *C‘, '■•=;•  14  1^:  • 
/  CAT  V  “  Y  /  L  C  2  4  3  f  :  l  ,  C :  T  4  2  ”  0  I 
7  c  3  ■* 

f  c  I  T  F  T  I  ■>  t  U  ~J  i  )  -  -  j '  .  1  = ;  j  ■  7  M  <  t  i_  i  .  ,  r 
/  FXCT  'nPT-(Yl2)»CPTi,'!:  =  L  ;  "  ' 

Cl  -i°  jt 

RH*«  =  10V  .  ,3C=-I.  ,11  =  1,  « 

IFL  =  3,L..ST  =  3, 

CE  NO 

f  CATV  "  r/,JG243".''2  »  *,C‘.  - 14  1  I 

/  CATV  VY/LC  24  3CI2,AC-!  ^.  '!4>f;  ; 

/  CATV  -  T  /  F  T  j  0  F  l  -  1  »  A  c  \  i  =  rTj>i  1 
f  t  2  J 


In  this  particular  case  year  5-7  are  input  to  the  preprocessor. 
Year  5  resides  as  the  fifth  file  on  #002022,  year  6  and  7  reside  as  the 
first  two  files  on  #000243. 


7 


Note  that  a  total  of  7  files  are  created  by  the  preprocessor  program. 
Each  year  has  2  files  (U  and  L)  and  FT30F001  contains  bottom  topography 
and  other  information. 

See  program  listing  for  interpretation  of  input  parameters. 

(ii)  Create  Load  Modules 


/  jcm  LI'l*G'T‘nL“  V)1  ,  1M  4‘jX.  ,L! ’T=(C,  ),  >) 
/  LI«1T  »4ND  =  5C  ,MlN=* 

/  PD  '!Y»USmCAr/077/JS. /L1..SLI 
/  A SG  *fl ,  •1Y/SC7R2lR1,US1  =VP. 

/  LNK  LNK0l»T  =  CN»r,A,S,F,.<l,.O,LSPACz*2i:?- 

LIBRARY  ’ll 
INCLUDE  ••nuULI 

/  CATV  N  Y  /  L  fl  .1 f)  1  »  A  C '  J  “  =  S  Y  G  «  L  !  M  ) 

/  EDJ 


/  J3B  LlNiGrULNND2,17A  4  S  A  ■  ,  L  1  N  r-L  i  , **P  T  =CG  . !) ,  O 
/  LIMIT  !3  A,ND  =  5G  =  : 

/  PD  MY  »USEP.C  AT /D77/«SC  /L.1  ">  JL  1 

/  ASG  D1 »  MY  /SCDR2LRD »USf  =  S.ir 

f  ASG  "Uli  J  »  M  Y /  NUR D ATEK  *  US  Z  =  S H 

/  LNS  LN!SMPT=(M,Y  .,/0,LSP  ',Ci  =  Z:r 

LIBRARY  n  I 

INCLUDE  KDDUL  2 

LIBRARY  .1  'i  J 

f  CATV  MY/LMTD2 »ACN *-S!  .!  ’  " 

/  EflJ 


The  small-core  version  of  the  2-layer  energetics  program  has  two 
driving  programs  (PROGRAM  MODULI  and  PROGRAM  M0DUL2).  These  two  programs 
should  be  executed  in  series.  The  above  two  jobs  create  load  modules  from 
the  object  library  and  store  the  modules  on  disk,  ready  for  execution  with 
appropriate  input  data. 
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(iii)  Execution  of  the  Load  Modules 


/  JOG  HG.G',1G.e:GZLFXUT,17A'4‘.A  ,  l  J  » l*.  i  ,  "»P  T  =C  C  ,  D  ,  R)  ,  OA  T 
f  LPIT  3ANO=5P 

/  PO  MY ,USfcRCAT/j77/^5./uI  l',Ll 
/  CD  CT  6  lPOul  * !'  A'iD  =  0/ 15/1  CP  ’*=  V  ;  ? 

/  FO  fttof-ooi  ,gagl>=i/e/:  ,rC- :=v::s 

/  FO  FT71FC0 lt'J'\.iO=l/^/».».''CFvi  =  VLj 
/  CD  =T72CC01»:H.'!0-1/Jfi»  <Cc  *  =  V'  'j 
/  A  SG  F  T  4  1 F  0  0  1 » **  Y  /  U  2  -  2  2  :G3,'JS  £=£.:" 

7  A  SG  F  T  4  2  F  0  0 1 »  0  Y  /  L  2  ■  2  ^  P  F  »  I S  E  =  S  i " 

/  A  SG  F  f  3  3  F  0  C  1  *  1 Y  X  "  T  3  2  .  j  1  *  !J  S  r  =  S  ir 
/  A  SG  LORD1  ,  n/L'^D,  ,'JS  =S«? 

/  FXQT  P  T  =  (  Y  » Z  )  » C P T  :  »  0  ^  =  L  C  0 1 

LMPUT1 

iyc=i »xsk:p=5,l*:ax=31  ,  ’33=:* ,t4$r  =  j,.:Ax-':LC  =  ii'  , 

IS  TH  =  1  ,  I'iUh  =  7S,.JST"  =  3»  .!T):i  =  vj,  I  J  *  <  -  1  r>  ^  ,  > 

CENG 

/  REL  FT41F301 
/  REL  F  T 4 2P  0  ■'  1 

/  ASG  F  T 4 1 1: 0 0  2  » *' Y / U  v 2 4  3 '  o  1 »  J S _  =  S  t  ^ 

/  ASG  FT42F001,  1Y/L2243  :  1*'JS:=S3' 

/  fkot  gpt=(y,Z),gpti»‘.e=9>*  >o.  ,g^=li';): 

GINP'JTl 

IYC  =  2  iiSKiP  =  0,L.‘  AX  =  36f  '3S=31  ,L  AST-3»  3  AX’,  CC  =  1  .  ^  . 
ISTH-l,I‘iDH  =  75»JSTH=3»  J  i.>’!  =  4  “■  j  I  J  A  Y=  i  5  '  / 

CEUO 

/  RFL  FT41F0C2 
/  REL  F  T 4  2 F 0 0  1 

/  ASG  FT41F(Jdl,MY/U0-,43-,:2*iJS£  =  S  P 
/  ASG  FT42F0D1,  1Y/LC243' 32,  JS:  =  S  l' 

/  FXOT  40T  =  (Y,Z),CPTi'!L  =  -'::  30 ,CP=l  ;<yo: 

SlrjPUTi 

IYC  =  3,ISKIP  =  0  ,Lr:aX  =  33,"  3i»*67  ,L  A 3 T  =3pP.  <f< EC  =2  3 
I S  T  H  =  1 ,  I  M')d=75  »  JST  !1  =  3  » 0.2  J'1  =  4  T  )  ’  f-l  5  .V  , 

CE  90 

/  P. c  L  FT4ir-0”i 
/  RFL  FT42F0n 
/  PEL  LOGOI 

/  ASG  L‘*3D2»MY/L0GD2,'JS-=S.l> 

/  ASG  FT41FC01 ,  M.r/U2022  *05»'JS0*S  jp 
/  ASG  FT42F001,YY/L2222  :2bjUSC  =  SrP 
/  FXGT  0  P  T  = ( Y  *  Z ) »  C  P  T 1 M  L - 9  0 j  •  J  v  »  G  G  =  L  P'  0  2 
6 1  N°  ’JT2 

IYC=1*  IS<IP  =  S,LKAX  =  j!  ,  -mXUEC  =  1:  :.,LASI  =  3,  IPr=?5»  , 

ISTH  =  1  ,  £'IUH  =  7S*  JSTM-3,  I  Ui.  =  4  », 

IDAY  =  150  .i  , 

GEflG 
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/  REL  FI41FLC1 
/  RtL  FT42F0  '1 

/  A S G  FI41FU'.'1,KY/!JJ243  Jl,  JS.r  =  S 
/  ASG  FT42F001  ,  *U/L0243  1,USE  =  S.I' 

/  FXOT  TPT  =  (Y,ZJ,CPilM;L=K*,J  j:  ,C."’  =  L  l',0: 
tlNP’JT? 

I Y  C  =  2  *  I S  <  1 P  =  C  » L  ‘■l  A  X  =  3  6  »  :  a  X  r,  E  C  =  1 G  j  *  L  A  S  T  =  3  *  I P  c  =  3  5  -  , » 
IS  TH  =  1 »  nDH  =  75,  JST.‘1=3,  J';.!J‘!  =  4o, 

IDAY  =  15(V.i, 

LENT 

/  REL  r T4 1 “0 F 1 
/  REL  FT42F0C1 

/  ASG  FT41F001  ,  "Y/U0243 ■  j-t  'JSE  =  S  iR 
/  ASG  FT42F001  * 'IY/LC'243  ■:  2, 'JS£=S  1r 
/  FD  FT11F0C1 ,3A‘.'0  =  4/2G/2 

/  FXOT  ftp  T  =  (  Y  »Z  ) ,  CPT1  ‘I t  =9  j 0  J »GC;  =L  1002  »  '.0D.’’2  "  =  24X 
cinpotz 

IYC  =  3.ISXIP  =  0,L  I  A  X  =  3  3  »  :AXKEC*ll  »  L  *.  S  T  -  3,IP  -=2f>»  j 
ISTH=l»I‘jOM  =  75|JSTM  =  3j  J.<3u=4?, 

ID  AY  =1S0  j  j 
CE'IO 

/REL  FT  41Ffl  1 
/  RFL  F  T  4  2  F  3 i 

/  cat  my/crass,  AC:n  =  pr:i.- i 

/  cftj 


This  is  the  main  program  of  the  two-layer  energetics  package.  The 
output  of  the  preprocessor  should  be  examined  by  the  user  to  determine 
the  time  interval  for  the  averaging  process.  For  this  case  the  first  five 
records  on  the  first  file  are  skipped  and  the  averaging  period  covers 
100  records.  See  program  listing  for  meanings  of  input  parameters. 
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(iv)  Switching  of  Storage  Device  for  Data  Files 


/  J  3 LlN.GDM.USMET.'U-L,  :/.<■  4fX.  ,LT  il-Li  t  )PT  =  Ct, 0,R),C  .T 
/  limit  <;a\p=5 

/  PO  MY,USfcP.CAT/!)77/iS  '/LI  .t'_* 

/  A  SG  «=l,  MY/U2G22c:5  ,'JS:  -S  :iT 
/  ASG  F2,mY/L2322F:.5>US  =  ST' 

/  ASG  r3,MY/UG243c:i»US  =  S.K 
/  ASG  F4, MY/L0243FG1 ,US  =S:;  l 
/  ASG  FS.MY/UC243F j2>US  =S  LG 
/  ASG  »  MY/L0243nu  2  »'JS:  =S.1U 
/  f  1 F ?,  FILE 

7  CATV  MY/'J2022c:'5,AC’M  =  f:,GTYF  =  T-,1>: 

/  CATV  f4.Y  /L202  2FG5  *  AC*I1  =  r  2  ,  IT  Y,1  =  T  f  *f 
/  CATV  .MY/IJC243F01  ,ACMP.  -F  2  »  D  T  YP  =  T  APL 

/  CATV  v  Y  /LG  24  3Fv  1  »  t  MM-P4  »  )T  YP  =  T  '  t 
/  CATV  MY/U0243FC  2,ACM  1=PS,DT  YP=T/:>f 
/  CATV  M  Y  /  L  <?  2  4  3  c '  2  »  A  C  lv.  -  r  0  » ; )  T  Y  P  =  T  A  ^  r 
/  MFRr- 
/  E^J 


/  Jflb  LIN.G1  *..TAPL2,1S'.,'.7>.',4  5\.  ,LI  m»l:,*>PT  =  CG,3,2),C'T 
/  L  Irt  I T  3Ar:0  =  42i.'»'j:  <  =  1 
/  PU  M  Y  » U  S  £  °.  C  AT  /  3  ?  7  /  3  S  ‘  /  L I  N  :•  L I 
/  PQ  F 1 , 3  A f , D  =  4 :) / 7  v / 5  ,  F  :> . •  vi  =  P S  ,  J .<  S 2  =  4 C  D  : ,  p  C F  •', *  V 
/  PD  c2  »!iA!il)  =  4P/7  0/  5  S2  =  41  ?  »  r.Cc"  =  ' 

/  FI)  F3f  f.AND  =  40/7G/5,F  )  :  j  =P  S  ,  r.K  S2  =4'  I  .  ,PCF”  =  V 
/  FD  r  4  »  2  A  NO  =  4  U  /  7  0  /  5  »  F  0  ->  G  =P  $  »  3 1*  S  2  =  4  C  j  C  »  2  C  P  "■  =  7  , 

/  FD  =5 ,3  VtO=4  V7P/5  ,-0  G  =  PS »  /..<S  I  =4f-T  •-'  *  ECr  !*V  v 

/  FD  F6 , 3  AND  *4P  /7 *V 5  , r  3-'  G  =l» S  ,  .U S 2  =40  ?  ,°.CFM  =  ,*iS 

t  MFR  FILE OUT 

/  ASG  FI  *  MY /U2C22FC  5  »US  =3  12 

i  /  ASG  F 2  »  1 Y / L 2 0 2 2 c C  5  » t J S '  =  S  i M 

/  ASG  F3,MY/UC243F01,,JS  =  b.L< 

/  ASG  F4»MY/LC243F01»US’  =  S,i'. 

/  ASG  FS,  3Y/UC243  "G2*US’  =S:iR 
/  ASG  F6,MY/L0243CC2,US.=S.1R 
/  HFP  = 

I  /  CATV  v  Y /U2 f  2 2 c(  5  »  A  C  V,  = 1 

/  CATV  MY/L2322FC  5  ,  ACGM-.-^ 

/  CATV  MY /UG243Fr  1,  Uh  i  =  F  3 
/  CATV  GY /Lo243cG  1 «  '.C‘i  1  =  .*4 
/  CATV  MY/U0243^J2» ACNM=Fb 
/  CATV  MY/LP243F'  2»  TC'-I- 

I  /  EDJ 


1 


11 


These  two  programs  switch  the  storage  device  for  the  data  files. 
To  avoid  high  disk  changes,  it  is  advisable  to  move  data  from  disk  to 
tape,  and  then  before  another  execution  of  the  energetics  program, 
move  them  back  to  disk.  Note  that  the  outputs  of  the  preprocessor 
reside  on  disk,  they  should  not  be  left  on  the  disk  overnight. 
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(v)  Program  Modification 


/  J11  L  IN,  ti  SCiJ'’  5M‘.»:  i>  .4r  »•  *  l  I  '*  t »_  i  *  *’  •’  1  =  l 

r  l  iMi  r  *  \  v»  =  ‘, '  ,  ‘i  •!  =  • 

/  Pi)  sr ,  jsr  rcat  /i>7  7/'  s r /l  •  .  li. : 

/  A  SG  S  3 ,  Rf  /SC  Lr.  S  .  VS  =  S  i  ' 

/  A  SG  "  i  ,  ,,r/SC)-!.''LR -t 'JS  =>.•' 

/  ClFr/ 

<<  cf'at  s  < ,sj 
<<  C^Pf  'M,'!.’ 

<<  f-S’.G::  SP  ,  T  7  Pf' ,  >)P  1)  A  1 

c  stc-cr  S'jrc 

-2S  7  » 2  6  S 

-2R  6 

IPCLPK3)  •  r:  J.  *>  C.LL  “ -  ~  I  F  C  i  ,  '  (  1 2  <•  ,  G  ?  "--"’V  )MTJ 

Rt;G'iC7r')  :: 

I  r  (  I  P  F  (  7  )  .  J  U  '  L P  L  "  T  *■  C '  1  »  *  (  v.  I  J  »  S  J  •’-.■I  •'>■•.  )  *  ,  1  M 

p  t  a  )  c  7  j )  . : 

I .  -  (  L  0  '  C  1 2  )  •  L.  <  »  .)  G  *  L ■•w"  r\  t*  :  ,  'CnP-  »c.*,  )'»!’> 

p-A  )(7G) 

!r(TPF(i7).c .)  C  '-’-  p  L  '■  r  l  C  -  !  »'(r>2LiG:  >  “' »  I  7  -> 

H  t  A  !  (  7  r>  )  ■!  1 

lr(  IPrC2i).  >  v  <L‘_  'J  -  -  r  l'  C  -  J  ,  *  C  <J  M  ,  •*  '  >  M  7  •> 

/  FT‘(  F  I  V  FKS  =  \X  *i  1  <"VT  =C  J  ,G  *  *  V  »  J  )  »  i  1-1 :  " 

/  f o  ft i6-;-  i ,  jr-i')  =  :’/i  /. 

/  CIFLR  C  A D -  R  =  i  ?  11 »  S t  u  - _  ^ 

<<  S°L  1 1  T i;  ■' "* »  S P 
<<  SPLIT  SYS.'V^^jG; 

<<  ppi  a  s? 

/  CATV  M /SCORJLRS, ACS  , 

/  CATV  “Y/SC'IRJIP  T»  ACS  "  -  : . 

/  f 


User  may  make  permanent  changes  to  the  package  by  executing  this 
program.  The  changes  should  be  bracketed  by  <<  MERGE  and  /FTN. 
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III.  IRREGULAR  GEOMETRY  CAPABILITY 


The  hydrodynamical  ocean  models  developed  at  NORDA  were  originally 
designed  for  rectangular  basins,  because  of  the  ease  in  applying  the  fast 
Helmholtz  solver  for  regular  geometries  and  the  simplicity  in  coding  the 
computer  programs.  Most  dynamic  feacures  of  ocean  circulation  can  be 
investigated  with  these  "square  oceans"  and,  indeed,  the  computation 
cost  is  relatively  low  because  of  the  fast  solver  and  the  high  degree  of 
vectorization.  However,  as  ocean  modeling  outgrows  its  infant  state, 
a  realistic  coastline  necessarily  becomes  a  major  factor  in  defining  the 
domain  of  the  flow  field. 

Historically,  the  successive  over  relaxation  method  (SOR)  is  used 
to  solve  the  Poisson  (or  Helmholtz)  equation  in  a  basin  of  realistic 
shape.  This  technique  is  time  consuming  (especially  when  a  large  number 
of  grid  points  is  concerned)  its  accuracy  is  not  quite  desirable,  and  the 
procedure  is  usually  difficult  to  vectorize.  In  order  to  extend  the 
modeling  capability  at  NORDA,  JAYCOR  proposed  to  implement  the  vectorized 
irregular  geometry  Helmholtz  solver  to  NORDA' s  existing  hydro-models. 

The  technique,  the  so-called  capacity  matrix  method,  was  developed  by 
Dr.  D.  R.  Moore  of  Cambridge  University.  It  has  been  successfully  used 
in  Dr.  H.  E.  Hurlburt's  free  surface  ocean  models  in  the  study  of  the 
Gulf  of  Mexico. 

The  implementation  is  done  in  the  following  way.  For  each  model 
we  have  two  packages  of  routines  instead  of  just  one  as  in  the  square 
basin  version.  The  first  package  is  almost  the  same  as  its  predecessor 
except  one  parameter  NREG  is  added  to  the  input  NAMELIST;  it  controls 
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the  basin  geometry.  If  NREG  equals  zero,  the  model  is  run  just  as  it 
would  with  the  regular  geometry  version.  If  NREG  is  not  zero,  it  stands 
for  the  number  of  land  blocks  added  to  an  otherwise  rectangular  basin, 
e.g.  NREG=3  means  3  land  blocks  are  present  in  the  basin.  When  NREG  is 
not  equal  to  zero,  a  number  of  subroutines  are  called  by  the  model 
program  to  activate  the  irregular  geometry  solver.  These  routines  reside 
in  the  second  package,  usually  named  with  abbreviations  of  Land  and  Sea 
such  as  "QGLNDSGS". 

The  present  irregular  Helmholtz  solver  can  handle  land  blocks  of 
arbitrary  size  (rectangular  though),  at  arbitrary  location  as  long  as  they 
don't  share  a  common  boundary  with  each  other.  Once  the  basin  geometry 
is  decided,  one  can  modify  the  subroutine  "IRGBDY"  in  the  second  package 
to  define  the  basin  shape  and  change  the  parameter  IBP  (for  number  of 
interior  boundary  points)  in  other  routines  of  the  same  package.  The 
first  package  need  not  be  changed  except  for  parameters  which  determine 
the  number  of  grid  points  (i.e.  ISF  and  JSF).  Followinq  is  a  detailed 
description  of  the  usage  of  four  hydro  models. 


A. 


Quasi -Geostrophi c  One-Layer  Model 
( i )  Run  Main  Program 


/  J*>i  Lit..  M  1  I  R  b-'  ,  :  7  X  '4S^.  »L  1  .  i  L  1  ,  =  T,  )  ,  C  .  T  =  1 

/  l  I  w  1  T  =  S  I 

/  PO  1Y  ,  IJStV.C  AT/ 17  7/ •  /l  • iL  I 

/  A  SG  :.Ll  »  11 S  .<  C  T  /  ‘ )  7  7/1  ‘j  /LI  - ,l.‘  ir=S  ! 

/  ASA  "1  ,  "Y  /  .  ,l!  Tr  “1  ,US  =  S  : 

/  A  SG  V,‘n/  lAASLVCl  ,  H\,;  =  S  I 
/  A  SC.  1.1  , -’.Y /'JOL*.  >SEO,'lj  =S.*T 
/  LV.  Lr.  •|PT  =  (i,  V,  A,S»r  »  .,<) 

LIUAr'Y  'll 
l  NCL'J  )£ 

LMRAKY  M 
LIBRARY  l 
L  1 6  R  A  :  Y  'll.  I 

/  FX'JT  •!PT=CY,:),CPT11‘;  --  t, ;  . 

GIN0  ) T 

XL  -  1  CO:  •  |YL  =  luCS  ».ClT  =  l/»i|  ,  .  ,  ,  A  -  3  .  •  If  !  T  =  S  >  .  ,  ‘ /  =  r  .  4  A  -  -  :  S  ,  „  =  1  .  b  <  i  -  1 1  . 

TN  V  =  -'.  S?.-r  A,  ;  =  ?.  f,2,  !  ' -JT  =  -  .  ,  P  ;  •  J  I-'i  .  ,  1  L  ’JL  '  ■: a  ,  I  x  =  A  6C  . , 

MAX-  )L  =  P  ,  >L  A  1  =  .’ .  S  ,  (  ,T.  :  1  -  -  L  .  » 

G  =  0  •  j  2  t  ;>  P  =  3  •  J  y 

V.  L  0  =  , , 

,‘R  E  G  =  1  » 

GE  NM 
/  c  TJ 


"QGL1TEMO"  is  the  updated  version  of  Lin's  QG  one-layer  model. 

A  parameter  NREG  is  added  to  the  input  namelist;  it  controls  the  basin 
geometry.  "QGLNDSEO"  contains  routines  bridging  the  irregular  geometry 
capability  to  the  Helmholtz  solver  for  regular  geometry  ("DANSLVRO") . 
See  listing  for  meanings  of  input  parameters. 
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(ii)  Update  Irregular  Geometry  Subroutines 


/  JO  3  LM.KX.OOUKEGM,:  7  O  .  4 0  a.  .  *  L  I  i  ,  'P  I  =  C ,  0 ,  '  ) ,  O  A  T  ~l 

/  L  r  «  I  T  «AN')  =  5r.  = 

/  PD  AT  ,  lj  S E  RC  AT  /  M 7 7-/;*  S  0  /  L I  'i  »L  1 
t  /  ASG  Sl»if/W»L’WSLS,!JS:  =  S :i  : 

/  ASG  01  ,  MY/  IGLNDSEO.’JS  =  S.H. 

/  CIFEK 

<<  CO^Y  SI  ,S2 

<<  COPY  31,02 

<<  MERGE  S2,TEil,'Jf>D;  Tl 

*  C  SELECT  CVKLtiJ 

-13,13 

COiiON/PfcRI/  IUH,  JoJY  ,?:>(»;•  U  \ 

-5C ,73 
-74,7A 

L'jTRT  I '  i  T  G  • ;  A  (X,l  ,  J  • ,  CS  T  ,  L  C'J  ,  I C  '«  .  ) 

»  REAL»R  ZOU ' 

-82,35 

60  ZDIJ  l=ZUUi+.<lC,Jj 

ZDU"  =  ZUUJ,+  "EPICS' 

Ii-Ctu  .  EU.  \  )  ;>.  T  j.<  , 

SIu  !•>=-/ J '  J  ‘ /  »Kfc  • 

'  00  70  J=i,J" 

no  ?o  i=i*i' 

70  X(I,J)=XCI,J)+S1G  v 

ZlR  O  =  0  . 

CALL  LA'lDCl  C  X  ,  1 .1 ,  i  0 ,  J  ! ,  Z  E  P r  ) 

<  SELECT  IkGoDY 

1  -1 

PARAMETER  I S  c  =  5  1 ,  i  S  c  =  '*  > 

-4 

CO  lO'j/PERI/  I :> .}X  ,  J-DY  ,  ^  r  *  Mi 

-10,10' 

Data  IC3/0, 1,1,0/ 

*  -11,11 

f>EK  =  C(J  i:Cl  )-2  )  +  C,:S;:-2  )*C  IS "-2  )♦(  TSr—  I  ,  »C  l)-D* 

1  <  I  MAC  U-2  >♦(  JSE--J  “I  CD- 17)  ’■:  ,5*2. 

ARE  A=CJMU1)-2J*C  TSr— 2)  +  (  ISE-I  :  \(i>-U*C  JSF-J1ICD3+P 
C  PICK  A  o.TUT  ON  r.-:c  T  VIE  LAM-Sc*  v^UiJARY 

IJUX=ISF-1 

‘  J3DY  = JSF 

<  SELECT  LAN9C 
-4 

1  ,  10  0(4, MAE Z) 

-14 

•w. 


I 
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entry  la.noci  (»  » *.  -a.  ,  .  r ,  ;r ,  vr  j'.  > 

D**  2  c  K=l,N2tG 
I  ■'  I  N  =  !  1 1  (  K  ) 

iKiti:$(A,K.)  .-.l.  )  : i  «  - i  1 1  co  + 1 

I  *!  A  X  =  I  *1 A  (  X  ) 

I P  (  I  tt '}  (  2  »  K  )  ,:u,  )  1HX  =  1 

jri*i=.  jncx  > 

it-(ib;U3,<)  .nl.  )  j  !K<)  + i 

J  1AX  =  J". A(f'  ) 

Ir(  IbHCl  ,*.)  •  '1l  »  )  J,t'U  =  J  If  U  )-: 

03  2t>  j=j'U*;.  j 
dt  2'  i  =  i^ri.  i  :,.x 
2:  WCI  ,  J)=VAL'Jr 

R  !;  T  >R  -J 

t  FIN  F  [  V  Er.  S  = f  i  X  »  F  T  :i  TP  T=C  J  ,A,  ,,Y  ,7  ,  j)  ,  i\  =  TF.  '  1 
/  C  I p  fc  P 

<<  SPLIT  TLMI,S.? 

<<  SPLIT  SY  S.V!''D,  ■}? 

<<  PRINT  S2 

/  CATV  »r/tJr,L  »JS2S»Wi  ‘  =  S. 

/  CATV  “  Y  /  ,'  G  L  N  )  S 1  ^  »  A  C  *1  ■•  =  J 

/  E  ^  J 


The  file  "QGLNDSGS"  (QG-LAND-SEA-SOURCE)  contains  the  irregular 
geometry  routines.  User  should  modify  this  file  to  suit  his  basin 
geometry.  Subroutine  "IRGBDY"  is  the  one  describing  the  basin  shape. 
The  above  is  for  permanent  changes  of  the  file.  All  changes  should 
be  bracketed  by  <<  MERGE  and  /FTN. 


I 
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I 


(iii)  Update  QG  One-layer  Model 


/  J03  H.-.IJXJASX  •  ,L  l-it'.i  ,-«‘T=(’w,  r,:  ),:M  '1 
/  LIMIT  )  =  5C  ,  1 1  N  =  _> 

/  po  my  ,'jst?.CAT/,'i77/:;b  j/l  i  .sl: 

/  A  SG  SI, MY /OGLITL  tS,’JS.  =S.!' 

/  ASG  01  ,  MY/'JCLl  TE  V>,'JS  .  =S.»1 

/  CIFEM 

<<  CO»Y  S 1  »  S  2 

<<  C’PY  71, M2 

<<  MFKGf  51  i  Tfc  ,'JPO^T  . 

<  SELECT  CALCU 
-13,13 

COMMON/PERi/  IuOAjJ  ,)Y 

-35 

il)U“.  =  O.UO'. 

IrC  MkEG  •  v  T •  C  l  'y'*  f  •>  4  2  3 
00  421  J  =  2,JSF>. 

1)0  421  1=2, 7 SF 
421  Zi)UM-=ZD'J-!*5(:,j,;.ll 

1°  £!1=C  ISr-2  )+CJS. —  2  )«■  1 
ZO'J'IsZOU'^IPERiSC  i  ,U>  ;1) 

SlGMA=-Z00'V(CiS.--2  X  C  JS!~-2  )  +  T  P  =  ) 
m  42  2  J  =  i  ,  J  S  F 
00  422  1=1, IS- 
42  2  SCI » J,  U)  =  SCI,J» 

GO  TO  424 

423  1 0  r<  E  =  1 

CALL  InT.jRA  CSCl,i,  1),IS.:,JS^,SC  iOX,J  .OY,  ;: 

424  CONTINUE 

<  SELECT  MAIN 
-98,93 

REAL*  9  ZOO" 

cwm/PERi/  ijOA.jijr 

-188 

1 3  D  X  =  1 

J  3  0  Y  =  1  2 
-237,237 

Z0'J3=U  *0  0" 

-241 ,242 

16  ZO!J  l  =  Zl)U  <  +  FS~.Cl ,  J  > 

Z:)UN  =  ZDUM*(ISP-2  )  +  C  JS--21+  1 
-245,246 

UERO^r 

CALL  INTuRA  ("S  ‘  ,  1  S  :  ,  JS1',"  ,7'>  »  ,  I  'is  "  ) 

-248,249 

RSa  =  l  ./Z  J(J‘. 


),2D'J  I"\E) 


) 
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v/  \/ 


-  33  J  *  3  3  , 

z  )'j*t=v..c.; . 

-334,336 

6  S  9  Z'J'J  i=Z  )U".-*-S(  ;  ,  J,  .<  ' 

z.)!M=zn!j  ’.*(c :5i  -z  )♦•( j  ',r-:  j  +  i  )•. : 1  ,  j  u , 

-33  7  *  33-' 

21  call  into-*  a  c  s  c  i » ' » o »  l  s  ~-  *  j  c  - ,  s  c :  ■  >  >. ,  j  .  i  y  , 

-340 , 341 

ZOU'!  =  Zi.HJ 
°.J.  :m  6  9'  *  Z  fvJ 
-343,343 

t  FT\  r  f  VE-sS  =  r.X  ,-l  *;-VT  =  C  i  ,  <■  ,  '  ,  Y  ,  ,  ))  ,  i  ,=T'.-  >» 

/  CIFE- 

<<  S  0  L I  T  TL9->,S' 

<  s° l i t  SYS..T'* r: 

<  PRIM  f> .. 


/  CATV  “Y/„,,L1T{-  !S,  L  :  -  0 

/  CATV  u.  Y  /GGL  1  T  ,  A  0 
/  E’J 


i. 


This  job  is  for  permanent  update  of  the  QG  one-layer  model.  All 
changes  should  be  bracketed  between  <<  MERGE  and  /FTN. 
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B. 


Rigid-lid  One-layer  Model 
(i )  Run  Main  Program 


"RGIINEWS"  is  the  updated  version  of  Lin's  Rigid-lid  one-layer 
model.  A  parameter  NREG  is  added  to  the  input  NAMELIST;  it  controls 
the  basin  geometry.  "LANDSEAS"  contains  routines  bridging  the  irregular 
geometry  capability  to  the  Helmholtz  solver  for  regular  geometry 
( "DANSLVRO" ) .  See  program  listing  for  meanings  of  input  parameters. 
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(ii)  Update  Rigid-lid  One-layer  Model 


.  JOB  LIN, 17X045X0, LIN$L1 ,OPT=(C,D,R) ,GAT=0 

/  LIMIT  dAMs5r>  ,  '1  l;=S 

/  PO  MY  »USERCAT/077  /u5i>/L1.iAL  i 

/  A SG  SI »  M  Y  /  RGL 1  *1E  w  S  » 'JS '"  =  Sm 

/  A  SG  31  ,  MY/RGLl:ifcW3,  US'  =S  I  3 

t  C  [  c  6  R 

<<  COPY  51, SJ 

<<  COPY  01,02 

<<  MERGE  S2  » T E Mp , UPJ AT L 

<  SELECT  CALCU 
-5,7 

-11,12 

COMMOU/NREG/  ;.UG 

-15,15 

1  0X2,  1 8  ,  T  R  D  X  2  ,  T  K !.) v  _ 

-43 

W 1  =  ■  j  . 

-70,7  > 

-76,83 

450  CALL  SOVLOHCM) 

-94,95 

-126,126 

503  CALL  LANDSC  C  I J  ,  G  .  ,  OX 2  , .-.D  Y 2  ) 

<  SELECT  MAI'. 

-94,97 

-101,101 

-106 

C0  riOM/AoC/AA,:i  i,CC 

-123,123 

1  0X2  , 1 3  ,  T R 0 X 2  ,  T «  T  c 

-177 

Tk 0X2  =  2.* 9  0X2 
TROY2=2  #*  RDYZ 
-225,22V 

IKMPEG  .GT.  GO  C  -LL  EV<L’.  ' 

/  ^  T  N  p  T  V  t  K  S  =  NX  »  F  T \ 9 ?  T=  (  J  , 6  ,  i ,  Y  ,  7  »  )),I  I  =  T  E  *  ° 
/  CIFEK 

<<  SPLIT  TEMP.S1 
«  SPLIT  SYS.OMOil,.:)? 

«  PRINT  52 

l  CATV  M  Y  /  R  GL  1 N  t  W  S  ,  A  C  N  *.  =  5  2 
/  CATV  M  Y  /  R  G  L  1  T  E  ■)  M  ,  A  C  M  •,  =  ! . 

r  e  o  j 


This  job  is  for  permanent  update  of  the  rigid-lid  one-layer  model, 
changes  should  be  bracketed  between « MERGE  and  /FTN. 
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(iii)  Update  Irregular  Geometry  Subroutines 


/  JflrJ  LI.N.1  7XC'.SXL  ,L1  Ni-Li  ,  H’  T  =  U  ,  T5  >F.  ) ,C  A  T  =1 
/  L  I P  1  T  -i AND  =  5 0  ,  = 

/  D  D  9  Y  ,  U  S L  RC  A T  /  3  7J  /  '  5  -•  /i_  [  ,  l  L  1 
/  A  SO  SI  ,  -1 Y  /  L  AN'T  S  t  AS  »  MS  =S  :  ■: 
f  ASO  91*  '1Y/LA:i')SfcAl,'jS  =  S  !  ’ 

/  C  I c  fc  R 
<<  C9PY  S  1 ,  S2 
<<  C'9PY  91,92 
<<  1ERGE  S 1 »  T E  ‘1 P  » I J P 0 A T L 
<  SELECT  I  i< r)  0  [)  V 
-10*10 

DATA  IJ3/U,l,:,r/ 

/  F  T  N  F  T  V  E  r,  S  =  NX  »  P  T  N  3  P  T  -  0  J  » .\  *  V  *  Y  ,'/»)),  I  N  =  T  E  P  P 
/  CI  =  ER 

<<  SPLIT  T  fc  MP , S  2 
<<  SPLIT  S  Y  S  •  3  R  9  3  *  0  2 
<<  PRINT  S? 

/  CATV  1Y/LA.<DSfcAS,4CV’=S' 

/  CATV  M  Y  /L  AMDS  R  A9  »  ACh  . 

/  E  9  J 


The  file  "LANDSEAS"  contains  the  irregular  geometry  subroutines. 
User  should  update  this  file  to  suit  his  basin  geometry.  Subroutine 
"IRGBDY"  is  the  one  describing  the  basin  shape.  The  above  job  is  for 
permanent  changes  of  the  file,  all  changes  should  be  bracketed  between 
<r<  MERGE  and  /FTN. 
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c. 


Quasi-Geostrophic  Two-Layer  Model 
(i)  Run  Main  Program 


i 


/  JOB  L  I  N  •  M  UL  •  G QL  2  N  E  h  i »  ?  7.\  '4bX  ,  L  I  1  ?  L  1  i'VT=CC#D»R)»C’«r  =  l 
/  LIMIT  GAG  1  =  5^  ,rtlN=v 
/  PD  MY  »  US  L  RC  AT/D  77  /  JS  j  /  L  I  i  *•  L  1 
/  A  SG  01  »  MY/UGL2Nt'.;9  ,'JS  '  =Si 
/  A  SG  92  »  M  Y  /U2L  DDSG’i  »  US '  =  S  I  ' 

/  ASG  0  3  »  ’  >  Y  /  D  A  S  L  V  P.  9  » U  S r  =  S M  ’ 

/  ASG  NL I ‘‘  t  !J S E R C A T /D 7 7  / L  t> i  / L I  *  U S c  =  S  i  . 

/  L NX  LNKnPT=C.,,Y,A,S,r,  j,>;) 

LIBRARY  91 
INCLUDE  '* AI  < 

LIBRARY  92 
LIBRARY  'M 
LIBRARY  fJLIT 

/  FXQT  9  P  T  =  (  Y  » Z  )  *  C  P  T  1  :  L  -  *  .•  j  j  . 

LINPUT 

XL  =  10P0.  ,YL  =  1P?  ;  .  ,  •)  r  =  144  ••  :u:=si  ♦  ,.=  y,.  .  ,GP  =  C.D2,U7',T  =  4 

FZ  =  7.43E-i  S,  .1  =  1 .  B.OE-11  , 1  ,''T  =  v  .  5_:_  ;,4  ,  l; ,  I  ,  M  TJ  T=  5  .  ,  *9  9D  T  =  5  2  .  ,  I 

TM  AX=36  i  •  »  M:\XCUL-2  •  AL  ■>  =  1  .  v  ,  \  3  6  =7  1 , 7  S  T  V.  T  =  , 

T  N  0  T  =  1  .  C  E-04,H‘JZ  =  »u  j  .  ,  '-P  =').,'  3  s  , 

TSTAKT=-l.,TB„X  =  6j.  ,  Y-:  j=  , 

P?.9UT  =  .K  .  , 

NR  EG  =  1 , 

LEND 
/  E  9  J 


"QGL2NEW0"  is  the  updated  version  of  Lin's  QG  two-layer  model.  A 
parameter  NREG  is  added  to  the  input  namelist;  it  controls  the  basins 
geometry.  "Q2LNDSGS"  contains  routines  bridging  the  irregular  geometry 
to  the  Helmholtz  solver  for  regular  geometry.  See  listing  for  meanings 
of  input  parameters. 
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in  f.j 


(ii)  Update  Model  Program 


/  J"H  LIU.  KDL.OGL<>  Ta.AGa.  »  u  1  » t  L 1  »  ^  P  T  =  t  C  » 'J  .f.  T 

/  LIMIT  n  V.">  =  5  ’  • 

/  PD  MV  , U$CKCAT/D77/'5j/L  !  \  i>  L  1 
/  cic::R 

/  a  so  si*‘'Y/aGL^‘:e<s,us  =$»’. 

/  A SG  31  ,'lY/:jGL2‘lPWO»Ub  =  S'i  : 

/  CIFEP 

<<  C^Y  S1.S2 

<<  CPpY  01,-12 

<<  MERGE  SI  ,  T E  •’ 0  ,  JPO  >%T  1- 


• 

/  FT**  F  T  V  E  K  S  =  f*X  »  c  T  U  0 !’  T-  (  U  » '. »  '*>  »<»  3)i  " 

/  CIFER 

<<  SPLIT  TciP.S? 

<<  SPLIT  SYS.  fl", -ID,  12 
<<  oRi-iT  S2 

/  CATV  M/ )GL2V1E'!S,«C v ‘  =SL 
/  CATV  MY/C'iL2NFW'1,  ALN  •-  •_ 

r  eoj 


This  job  is  for  permanent  update  of  the  QG  two-layer  program. 
All  changes  should  be  bracketed  by  <<MERGE  and  /FTN. 


t 


t 


i 
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I 


(iii)  Update  Irregular  Geometry  Subroutines 


i 


i 


/  JOB  LIU.  MUL.JOLPLElS'  ,  .  7  >.  ,  L  '  \  Hi  ,  T  =  (C  ,  0  ,  R  )  ,CA  T  n 

/  L  I u I T  dArjD=5>i'  ,  MlN  =  ‘. 

/  PD  MY  ,USERCAT/D77 /oSI /LI  ,'.LI 
/  A SG  Si  ,  MY/'JGLND SE S » US’  =S:i  ’ 

/  ASG  0 1  ,  M  Y  /  OGL  NO  SE  0  ,  :j  3 '  =  3  i1’. 

/  CICEP. 

<<  COPY  S 1 ,  S2 
<<  COPY  01,02 
<<  MERGE  SI  jTEMM  ,'JPDAT  - 


/  FTN  cl  VEkS  =  UX  ,FT  HPT  =CJ  ,  v  ,  Y  ,  V  ,  1)  ,  I  j  =  l t*  *' 
/  CIFGR 

<<  SPLIT  TERM, 32 
<<  S  °  L I  T  SYS. 'V!*!,'),  •’2 
<<  PRINT  S2 

/  CAT  MY  /J2LNDSES,  AC‘.  '.-G: 

/  CAT  MY/u2Ln;)SE0,  ;l-.  --  . 

/  E  0  J 


I 


This  job  updates  permanently  the  land  sea  package  for  a  particular 
basin  shape.  Replace  dots  with  your  changes. 


1 


I 


I 


26 
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D.  Rigid-lid  Two-Layer  Model 
(i)  Run  Main  Program 


/  JflD  Ll:..“OL.Rr>L?licW'J,i7X,ASx.  .LlMLi  ,  ^  P  T  =  t  C  ,  0  »  A  )  ,  C  ’*  1  =c 
/  LIMIT  rtAND=5D ,11N=1  ■' 

/  PD  »ir  »l)SrRCAT/r)77/D5L‘/Li  J»L1 
/  ASG  <H»«Y /RGl.VlEWOtUS!  =  S  I  . 
t  ASG  C^RY/RZLNDSE'D.'Jo  =  S‘!3 
/  ASG  03,  RY/DAiSLVRO  |MS'-=  S  i-' 

/  ASG  HH"»OS£‘>CAT/n77/lW_:  .»IJ3r=S;!  ’ 

/  ASG  S 1  »  3  Y  /  ,<  2  L  D  D  S  l  S  » 'J  S '  =  S  i 
/  CIF6P. 

<<  A  £  A  G  fc  Si  »7E'‘r’»,Ji,r'«T- 
<  SELECT  I.AG3DY 

-7,3 

DATA  J-'.I/l/ 

DATA  .m/15/ 

-10,1  » 

DAT  A  i3:J/l,l,i-,  / 

/  FTN  F  T  V  E  R  S  =  F  X  » E  T  4  ^  P  T  =  0  J  ».*  »  > » Y  ,7  ,  ) ) ,  T '  i  =  T  E  ‘  ° 

/  IP*  LMCT°T=(  ,,r,  A,S,  .  ,  O 
library  -j: 

INCLUDE  !t  M  ’< 

L  13  PART 
LIBRARY  Cl  3 
LIBRARY  DLLs 

/  FXOT  *1PT  =  (Y,Z),CPTI‘t£=y-.  V. 

LINP'IT 

XL=  1  .  ,  •  Y  L  - 1 0  0 '  )T  =  ic  =  \  =  3  D  J . , G  P  =  '  .  :  2  ,  MT  V  =  A  n  :) .  , 

FZ=7.4  3.  -'  5 , 3=1.38£-11  ,T:;U  =  j.5E-  ,G=?.8»  1..  =  :  ,.l;  7  JT  =  5  .  ,PR.VIT  =  9; .  ,  I£ULtR=43D, 

T  3  A  X  =  3  6  .  ,,!AX£DL=2,Am  =1.' ,  T  ST  A-?  T  =  o.  , 

or =120 » ni)2=dLM>.  ,gp=:  .  ;>  35  ,T.M-tT  =  1  r-  4 , 

TMAX=3:  .  ,TSTAD  T=-l  .  ,  lA  j>=l*PRM'T  =  i£«, 

T M  A X - 6 1  •  •  ,  PR  O’JT  =  3  3  •  , 

A=ioa: . , 
tE  MO 
/  EOJ 


This  job  run  is  a  rigid-lid,  two-layer  model  with  a  flat  bottom 

> 

and  closed  basin.  "RGL2NEW0"  contains  the  main  calculation.  "R2LNDSE0" 
contains  land-sea  information. 

» 


I 
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(ii)  Update  Irregular  Geometry  Subroutines 


/  JrtK  LI*UKOL.RGL2UP,_U,17X' ,L1  ,v_l#'i,'T‘=CC»'),R)#C*.  T  =  1 
/  L  I M I T  UA.V^Sj  = 

/  P  D  MY  » 1)  S  t  R  C  A  T  /  U  7  7  /  b  5  C  /  L !  N  ?•  L  1 
/  A  SC.  SI » 1 Y  /  R2LND  SE  S  » ‘JS  ;  =  SMT 
/  ASS  01  »MY/R2L'.'!)SEn,'JSt  =SM;'. 

/  CIFER 

<<  COPY  S 1 »  S  2 

<<  COPY  01, '12 

<<  MERGE  S2 »TE  1°, UPDATE 


/  FTN  FTVERS  =  NX,r-T!l';rT  =  <lJ,N,M,Y  ,V,  J),  :  ;  =  TEl  1> 
/  CIFER 

<<  S°LI  T  TEM2,S2 

<<  split  sYs.TnD,o: 

<<  PRIM  S 2 

/  CATV  MY/ H2L NOSES, AC'  * " S . 

/  CATV  MY/R2LH'3SEM,ACN  i-  i 
/  ECJ 


This  job  is  for  permanent  update  of  the  land  sea  package. 
Changes  are  inserted  at  the  dots. 
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(iii)  Update  Model  Program 


/  J  '*B  L  1  ■  • .  'i  D  L  .  R  G  LZ  U  n  J  T  t  IT  •  4  X  .  _  t  ( *  L  1  ,  m '*  T  =  t  C  ,  .I  ,  )  ,  CA  T 
/  LinT  0  A-iO  =  5G  ,  riN=- 
/  PD  MY  »  USE  SCAT/. 177/u  5  .  /Li  ^  L  1 
/  A  SG  SI  »  MY /RGL2NE  ^  S  »  MS  c  =  SH ’ 

/  A  SG  01»MY/RGL.?ilEiJT»')S'=S  !’ 

/  CIFLR  C  \Or>MC“  =  J'.,K  ,S~CnZ 
<<  CMPY  S 1 »  Sc 
<<  C5°Y  Ml,*!2 
«  MERGE  SZ , TEMP , UPDATE 


/  FTN  F  T  V  E P,  S  =  NX  »  F  T  ’I ’TP  T  =  t >J  »  \  »  '• »  V  >  V  ,  j)»IN  =  lfc‘‘P 
f  CIctP.  CAUOiE  >  =  .i  ?  K  » S  E  C  =  i  2 
<<  SPLIT  TfcMP,$2 
<<  split  SYs.rriu,  j; 

<<  PRINT  S2 

/  CATY  NY  /  k  G  L  Z  '•  E  •»  S  t  A  C  •"!  1  -  S 
/  CATV  NY /«GL2,<E'H»;.L'i  =  ‘_ 

/  EMJ 


This  is  for  permanent  update  of  the  main  calculation  program. 
The  changes  should  replace  the  dots. 
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IV.  SEAMOUNT 


A  joint  study  is  presently  being  conducted  with  Dr.  Andrew  Vastano 
of  Texas  A&M  University  on  flow  effects  near  a  seamount.  Dr.  Vastano 
has  gathered  data  from  drogues  released  in  the  western  North  Pacific 
Ocean.  The  drogues  passed  through  the  Emperor  Seamount  chain  located  at 
approximately  40°N,  170°E.  The  path  of  the  drogues  was  dramatically 
altered  in  the  vincinity  of  the  seamounts  (see  Fig.  1).  Using  models 
presently  available  at  Code  322,  NORDA,  the  dynamic  effects  of  flow 
around  a  seamount  are  being  investigated. 

The  first  model  chosen  for  use  in  this  study  was  the  barotropic 
model.  This  model  demonstrates  the  effect  of  the  external  mode.  The 
model  is  designed  to  have  a  steady  prescribed  inflow  along  the  entire 
western  boundary  and  an  open  eastern  boundary  with  conservation  of 
mass  transport  strictly  enforced.  Numerous  experiments  were  made  using 
10  km  grid  spacing  in  x  and  y,  a  coefficient  of  eddy  viscosity  of  3x10^ 
dynes/cm2  and  an  inflow  velocity  of  10  cm/sec.  When  the  model  was  run 
with  no  topography,  solutions  were  steady  out  to  90  days.  However, 
when  a  small  magnitude  seamount  was  introduced,  a  disturbance  appeared 
on  the  upstream  side  of  the  seamount.  Tests  showed  that  the  steepness 
of  the  seamount  was  not  the  cause  of  the  instability.  The  disturbance 
disappeared  when  a  westward  flow  was  used  or  when  beta  was  set  equal  to 
zero.  The  western  boundary  condition  does  not  allow  the  "Rossby  Wave” 
out  of  that  boundary  and  the  solutions  degenerate  (Fig.  2). 

A  number  of  methods  were  tried  to  alleviate  this  problem  at  the 
western  boundary.  It  was  finally  decided  to  use  viscous  damping  at  the 
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boundary  of  the  form  Cg  V  where  is  a  drag  coefficient  of  5xl0‘‘ .  This 
boundary  damping  was  also  used  at  the  outflow  boundary  to  avoid  any 
instability  on  outflow.  Figure  3  shows  that  the  solution  using  the 
boundary  damping  ran  to  at  least  200  days  without  the  instabilities  of  the 
previous  cases. 

After  testing  the  effects  of  the  location  of  the  seamount  in  a 
basin  the  size  of  that  in  Fig.  3,  and  a  test  with  a  topography  like 
that  of  the  Emperor  Seamount,  a  two-layer  model  will  be  used  and  drifter 
tracks  will  be  made  from  those  solutions  to  compare  with  the  real  data. 
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REE  SURFACE  DEVIATION 
IH=  0.05  M  DAY=  40 
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